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REDUCED TWO-TYPE DECOMPOSABLE CRITICAL BRANCHING 
PROCESSES WITH POSSIBLY INFINITE VARIANCE 

CHARLINE SMADI AND VLADIMIR A. VATUTIN 


Abstract. We consider a Galton-Watson process Z(n) = (Zi(n), Z 2 (n)) with two types of 
particles. Particles of type 2 may produce offspring of both types while particles of type 1 may 
produce particles of their own type only. Let Zi (m, n) be the number of particles of type i 
at time m < n having offspring at time n. Assuming that the process is critical and that the 
variance of the offspring distribution may be infinite we describe the asymptotic behavior, as 
m, n —> oo of the law of Z(m, n) = (Zi(m, n), Z 2 (m, ri)) given Z(n) ^ 0. We find three different 
types of coexistence of particles of both types. Besides, we describe, in the three cases, the 
distributions of the birth time and the type of the most recent common ancestor of individuals 
alive at time n —>■ oo. 


1. Introduction 

Describing the genealogy of populations yields a better understanding of patterns of variation 
in genetic data. It is a key issue in population genetics and has been the subject of many 
investigations in recent years. In this work, we are interested in the genealogy of a population 
modeled by a two-type decomposable Galton-Watson branching process conditionally on its 
survival until a large time n. 

Let Zi(n) be the number of type i particles (i € {1,2}) alive at time n E {0,1,...}, Z (n) = 
(Zi(n), Z- 2 (n)), and denote by Z(0) = e ? ; the event 

Zi(o) + ^ 2 ( 0 ) = Zi( 0 ) = 1 . 

We assume that the process is critical, that is to say 

E[Z,(l)|Z(0)= ei ] = l, *€{1,2}. 

Type 1 particles may produce only particles of their own type while type 2 particles are able 
to produce particles of both types. The production of type 1 particles by type 2 particles can 
be seen as recurrent mutations or migrations. In the finite variance case, this process has been 
studied in nans], and in the case of infinite variance it has been investigated in [28j and i 2 11 
in detail. Zubkov [28] and Vatutin and Sagitov j2lj have analyzed the asymptotic behavior of 
the survival probability of the process given Z(0) = e 2 , and described the conditional limiting 
distribution of the population size at time n —>• 00 given Z(n) / 0 (where 0 = (0,0)) and 
Z(0) = e 2 - Our motivation for considering the infinite variance case comes from [ 8 ]. In this 
work Eldon and Wakeley, using genetic data from a population of pacific oysters, showed that 
one individual was able to generate enough offsprings in one birth event to produce a significant 
proportion of the population. This is also the case for many types of fungi and viruses. 

The present paper deals with the structure of the family tree of the process {Z (m), 0 < m < n} 
given Z (n) 7 ^ 0. More precisely, we investigate the properties of the process 

Z (m,n) = (Zi(m, n), ^(m, n)) 

where Z,- (m, n) is the number of type i particles alive at time m < n and having descendants 
at time n. The process Z (-,n) is called a reduced branching process and can be thought of as 
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the family tree relating the individuals alive at time n. Two important characteristics of the 
reduced process are the birth moment fi n and the type T n of the most recent common ancestor 
(MRCA) of all individuals alive at time n defined as 

fin = max {m < n : Z x (m,n) + Z 2 (m,ri) = 1 } , T n = {i G {1, 2} : Z* (fi n , n) = 1 }. ( 1 ) 


Limiting reduced trees have been extensively studied in the literature. In the monotype and 
finite variance case, the distribution of fi n already appears in Zubkov m, and the full structure 
of the reduced tree is derived by Fleischmann and Siegmund-Schultze [ 8 j. Analogous results 
in the stable case and/or in the more general setting of multitype indecomposable branching 
processes with finite or infinite offspring variance can be found in Vatutin [161 20] and Yakymiv 
[25]. Finally, Duquesne and Le Gall [5] derived general properties of limiting reduced trees in 
the continuous setting. 

For si, s 2 € [0,1], let 


Fi(s i) = E 




Z(0) = ei and F 2 i(si, s 2 ) = E 


zfi 1)Z 2 (1) 


Z(0) — e 2 


be the offspring generating functions for particles of both types. Our basic assumptions on the 
offspring distributions are formulated in terms of the generating functions F\ (sq) and E 2 i(si, S 2 ) 
as si, s 2 t 1 ; 

Fi(si) = si + (1 — si) 1+ai L\ (1 — si), (2) 

and 

F 2 i{si,s 2 ) = s 2 + (1 - s 2 ) 1+ " 2 L 2 (1 - s 2 ) - (A 2 1 - p(si,s 2 )) (1 - si), (3) 

where 0 < or, a 2 < 1, L\(x) and L 2 (x) are slowly varying functions as x | 0, 

A 21 = E[Z 1 (l)|Z(0) = e 2 ]€(0,oo), (4) 

and the function p(si,s 2 ) -A 0 as sis 2 f 1- The form of the probability generating function F\ 
is natural if we aim at modeling a population where the offspring distribution has an infinite 
variance. Indeed, ([2]) holds if and only if the distribution of Z\(l) is in the domain of attraction 
of a stable law with index 1 + or. Such a monotype branching process has been investigated by 
many authors (see for example ca m Eg Ea). 

Introduce the following notation, for i € {1, 2}: 

Qi(n) = P (Zi(ra) > 0|Z(0) = e*) and Q 2 i(n) = P (Z(n) / 0|Z(0) = e 2 ). 

It is known m that, given Condition ([3]) with si = 1 and Condition ([2]), 

Qi(n) = (n + l) - 1 /“Tj(n), n -A 00, ( 5 ) 

where li(n) is a slowly varying function as n A 00, 

( 6 ) 


and for A > 0, 


lim E 

n—>00 


ainQ^{n)Li (Qi(n)) ~ 1, n -A 00, 


g-XQfi^zfin) | Z (Q) = ei ,Zi(n) > ol =! — (! + 


(7) 


We will see that three types of behavior emerge for the reduced trees, depending on the 
relative values of ot\ and a 2 . More precisely, we will distinguish three cases: 


nQi(n) = o (Q 2 (n )), n -A 00, (Theorems Q] and 0]) , ( 8 ) 

Q 2 (n) = o(nQi(n)), n -A 00, (Theorems [2] and [5]) , (9) 

nQ\{n) ~ crQ 2 (n), <r G (0, 00), n ^ 00, (Theorems [3] and [ 6 ]) . (10) 

We will show that under Condition ([9]), the reduced processes have similar properties as in the 
case of finite variance studied in detail by Vatutin in [20] and [23], while under condition ([ 8 ]) 
or m the dynamics of the reduced processes has an essentially different nature. 
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The paper is organized as follows. In Section [2] we describe the local behavior of the reduced 
process. Section [3] is devoted to its global behavior and to the law of the MRCA of individuals 
alive at a large time n. In Section [4] we give auxiliary results needed in the proofs. Sections [5] 
to [7] are devoted to the analysis of the local structure of the reduced process. In Sections [8] and 
[9] we consider the global behavior of the reduced process. Finally, in Section fTOl we derive the 
MRCA limit distributions. 

In the sequel P, and E, will denote the probability and expectation conditionally on Z(0) = e*. 
Sometimes it will be convenient to write P and E for P 2 and E 2 . Finally, the relations 

a(n) ~ b(n), a(n ) = 0(b(n)), a(n) = o(b(n)), and a(n) <C b(n) 

will be understood (if otherwise is not stated) as n —» oo and the symbols C, C\ , C 2 , • • • • will 
denote positive constants which may vary from line to line. 


2 . Local characteristics of the limit processes 


We first deal with Condition (|8j). Notice that a\ is necessarily less than 1, as otherwise nQ\{ri) 
would be a slowly varying function while Q 2 ( 71 ) is a regularly varying function with negative 
index. This allows us to introduce an integer-valued function h*(n ) satisfying 


a 1-A21 
1 — 


/i*(n)Qi(/i*(n)) ~Q 2 (n). 


( 11 ) 


Using properties of regularly varying functions (see, for instance, El, Ch.l, Lemma 1.8), we get 

h*(n) = n ai/Ml ~ ai)) l h ,(n), (12) 

where ^*(n) is a function slowly varying at infinity. Moreover, h*(n ) <C n (see Lemma [J 
page [9]). In what follows equalities of the form Z (x,n) = y for x,y € (0,oo) will be understood 
as Z([x] , n) = [y] where [w] stands for the integer part of w. The same agreement will be in 
force in other similar situations. For instance, Qi(x),x € [0,oo) will be treated as Qi([x]). 


Theorem 1. Let Conditions 
(0) If 1 < m <C n, then 


lim E 

n —> 00 


P|) and © be satisfied and (s\,S 2 ) be in (0, l) 2 . 

Zi(m,n) Z2(m,n)i 


|Z(n) + 0 


= s 2 ; 


(1) For any a € (0,1) 


lim E 

n—>00 


Zi ( an , n ) Z2(an,n)i 


Z(n)/0 =1- ( a + (l-a)(l- S2 )-“ 2 ) “ 2 


(13) 


(14) 


(2) For any function h satisfying h*(n) <C h(n) <C n and m = n — h(n), and A 2 > 0 


lim E 

n—too 


^Zi(m,n) 


exp < —A 


Q 2 (n) 


Q 2 i(h(n)) 

(3) For any t > 0, A 2 > 0 and m = n — th*(n ) 

Q 2 (n) 


Z 2 (m,n) > |Z(n) / 0 


= i-(i + a 2 -“ 2 ) 


-a 2 \ _1 / a2 


lim E 

n—>00 


Zi ( m , n ) 


exp < - A; 


‘Q 2 i(xh*{n)Y 


Z 2 (m,n) > |Z(n) / 0 


= 1 — ( 1 + (V~ 1/ai ( 1-si 


. 1-01 


+ A 2 ^j 


—OL2 


— l/a2 
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(4) For any function h satisfying 1 <C h(n) <C h*{n), m = n — h{n), and Ai, A 2 > 0 


lim E 

n—>00 


, , Q\{h*(ri)) Q 2 {n) Wn 

exp < ~ Al wwr Zl(m -" ) ‘ A2 feM5) Z2(ro '" ) J Z(n) # 0 


-i (i + bV+C " S J 


— 1 / 0:2 


We see that under the conditions of Theorem Q] the reduced process essentially consists of one 
particle of the second type at the initial stage, it consists of many second type particles but still 
no particles of the first type at the intermediate stage, and particles of the first type appear in 
the process only at time m = n — th*(n), t > 0 . As a consequence, the coexistence of both types 
in the reduced process is possible only within a relatively small (in comparison with n) time 
interval of order h*(n), located at the end of the time frame under consideration, [0, n\. Thus, 
the evolution of this reduced process is mainly supported by particles of the initial type. 


Let us now deal with Condition ©. To formulate the desired result we need to introduce a 
function (f> = 4>(Xi, X 2 ), Ai, A 2 > 0 , which solves the partial differential equation 


(1 + a 2 )Xi—^- + A 2 WT- — 4> — <f 1+a2 + a 2 A 2 iXi 
uX 1 0 X 2 


(15) 


with initial conditions 


0 ( 0 , 0 ) = 0 , 


d(f> (Ai, A 2 ) 


I Ai=A 2=0 '1 


<9Ai |A1=A2=U '* 21 ’ dX 2 

The problem of existence and uniqueness of a solution for Equation (11511 has been solved in 
and m- Note that if or = 0:2 = 1, the solution of (1151) has an explicit form: 

, /, . \ r~i —r - A 2 + \M 21 A 1 tanh \M 21 A 1 

* < 11 2) = ' /A2lAl A 2 tanhVA 21 Ai + V^lA,' 

Let < 7 * (to) be an integer-valued function satisfying 

Q 2 (g*(n)) ~ Q 2 i(n). (16) 

Using properties of regularly varying functions and the asymptotic representation for Q 2 \(n) 
given later on in (|34l) one may check that 

g*(n) =n a2/ ( ai(1+a2)) l g *(n), (17) 

where l g *(n ) is a function slowly varying at infinity, and that g*(n ) <C n. 

Theorem 2. Let Conditions and m be satisfied and (si,S 2 ) be in ( 0 , l) 2 . 

(0) I/l<m«g*(n), then 


d(f> (Ai, A 2 ) 


|Ai=A 2=0 1 - 


lim E 

n—>• 00 


Zi(m,n) Z 2 (rn^n) 


Z(n) + 0 


= S 2 ; 


(1) For any t > 0 


lim E 

n—>00 


Z 1 (tg*(n),n) Z 2 (tg*(n),n) 


Z(n) ^ 0 ] = 1 - r 1 / 02 </>((! - si)^^, (1 - s 2 )t 1/a2 ); 


1 °2 
(2) If g*(n) <C m <C n, then 


02^21 


lim E 

n—> oo 


s Z l( m,n) s Z2(m,n)| Z ( n ) + Q 1 = j _ (1 _ ^l/fl+oa). 


(3) For any a € (0,1) 


lim E 

n—> 00 


sf lM S ^an,n)^ n) ^ (jl = X _ ( a + (1 _ o) (l _ Sl )-«l) -V(«l (l+a 2 )) 


(18) 
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(4) If m = n — h(ri) for some 1 <C h(n) <C n, then for Ai > 0 


lim E 

n—too 


exp 


—Ai 


Qi{n) 

Qi{h(n)) 


Zi(m, n ) > s', 


Z 2 (m,n) 


Z(n) + 0 


= i-(i + Ar“ i ) 


-1/(01(1+02)) 


Once again, coexistence of both types in the reduced process is possible only within a rela¬ 
tively small time interval (here of order g*{n)). But unlike the case of Theorem □ this interval 
is located at the beginning of the population evolution and the reduced process is essentially 
composed of type 1 individuals. 


The next theorem deals with Condition (fTOl) . Let h E (l,oo) be the unique positive solution 


of 


-1+Q2 


— X = < 7 a 2 A 2 l, 


(19) 


and if = if (Ai, A 2 ), Ai, A 2 > 0 , the unique solution of the quasi-linear partial differential equation 


dif 


dif 


(1 + a 2 ) Ai—— + X 2 —— — if — b a2 if 1+a2 + a 2 A 2 i\i ( 1 + ( —Ai 


<9Ai 


‘dX 


a 




with initial conditions 


/. cn rp _ n ^(Ai,A 2 ) _ dif (Ar, A 2 ) , _ 

Y (9j UJ U, -- |a 1= A 2 = 0 — A 21 , -7^- |Ai=A 2 =0 — !• 

Then we have the following asymptotics for the intermediate case (HOD : 

Theorem 3. Let Conditions and m be satisfied, and (si,S 2 ) be in (0,1) 2 . 

(0) If 1 < m <C n, then 


( 20 ) 

( 21 ) 


lim E 

n—>00 


Z\ ( m,n ) Z2 ( m,n ) i 
S 1 S 2 


Z(n) / 0 


= s 2 ; 


( 22 ) 


(1) For any a € (0,1), 


lim E 

n—>00 


Zi(an,n) Z2(an,n) 
3 1 S 2 


Z(n) / 0 


(2) If m = n — h{n), where 1 <C h(n) <C n, then for Ai, A 2 > 0 


l/a 2 \ 


lim E 

n—>00 


exp 


XiQ 2 i(n) . . A 2 Q 2 (n) . U i 7n , n 

-Zi(m,n) — Z 2 (m,n) > Z (re) 7 ^ 0 


— 1 — if (Ai, A 2 ) • 


nQi{h{n)) ’ Q 2 {h{n))‘ 

In this last case, both types of particles coexist during an interval of order n. 

As we will see in Section [3l Theorems □ED will be elements of the proofs of more general 
statements dealing with conditional functional limit theorems for reduced branching processes. 


3 . Global description of the limit processes 

In this section, we study in detail the dynamics of the reduced process, in particular, the 
transitions between monotype and two-type phases. Let x 1 ^ = x(x — 1 ) • • • (x — k + 1 ). In 
Definitions m we specify five auxiliary continuous time Markov processes needed to describe 
the transitions. 

Definition 1 . We denote by {v = {y\, v 2 ), (<71 (si), ii±), (521 (si, s 2 ), /12)} a homogeneous Markov 
branching process initiated at time t = 0 by particles of type 1 and v 2 particles of type 2 . At 
rate //1 (resp. g 2 ) a particle of type 1 ( resp. 2) dies and produces children in accordance with 
the probability generating function g\ {resp. g 2 \). The newborn particles behave as their parents 
and independently of each other and of the past. Notice that when type 1 {resp. 2 ) particles 
do not evolve (no death, no offspring), g\{sf) = si {resp. g 2 i(si,s 2 ) = s 2 ) and we can take 
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any g\ ( resp. /i 2 )■ We will write gi, i € {1,2} to make it clear in the definition. With this 
representation in hands we can now introduce the following four branching processes: 

( 1 ) {X(t),t > 0} = {e 2 , (si,/xi), (52 A) ('Si,S2), 1 )}, where 


g { 2 X \si, 82) = —((1- s 2 ) 1+ “ 2 - 1 + (1 + a 2 )s 2 ) = - f; (1+ ,^ 2)W (-1)M 
a 2 «2 k\ 

k=Z 

(2) {Y (t) = Y 2 {t)) ,t > 0} = {e 2 , (si,/xi), (g%i\si, s 2 ), 1 + l/a 2 )}, where 

g 2 i\si,S 2 ) = —Y((! - s 2 ) 1+ “ 2 - 1 + (1 + a 2 )s 2 + «i) 


+ a 2 

1 

1 + a 2 


I_V' i)k s k + 1 „ _ 

- k\ 2 1 + a 2 


k=2 


( 3 ) {V(t),t > 0} = {ei, (gf } (si), 1), (s 2 ,g 2 )}, where 


g { i\s 1 ) =--[(1 - si) 71 +7i«i - 1 

71 - 1 L 


1 00 [ fc ] 

Zfri-Yst 


71 — 1 ' /c! 

fc =2 


(23) 


(24) 


(25) 


and 71 = 1 + «i(l + a 2 ). 

(^ {W(f) = (Wi(t), W 2 (t)) ,t > 0} = {e 2 , (g ( ' l \si), 1), (g^ (si, s 2 ),«)}, where 

(1 + a 2 ) 6“ 2 - 1 


k : 


0.2 


> 1 , 


ffr ) (^i) = -((l-^i) 1+ “ 1 -l + (l + 

Oi\ 1 


1 y-(l + Ql) [t| 




Q7 


k='2 


k\ 


(26) 

(27) 


and 


(W), x 1 /'<ra 2J 4 2 i 

521 ( s i> s 2) = 


si + 6 Q2 ((1-s 2 ) 1+ “ 2 -1 + (1 + 


a 2 A V 

^i Sl + V ( 1 + « 2 ) [fcl (_i )fca fc. 

bn a 2 « 2 —' 


fc =2 


fc! 


(28) 


Note that X(-) and Y(-) are monotype processes, while Y(-) and W(-) are two-type processes. 
Moreover, particles of type 1 of the process Y(-) are sterile and immortal. 

Definition 2. We denote by |G(f) = (Gi(i), G 2 (t)),t > 0} a two-dimensional Markov process 
with values in NoX R + = {0,1,2,...} x (0, 00). The initial state of G(-) is random: 

G(0) = (G 1 (0),G 2 (0)) = (O,0 2 ) 

where the distribution of O2 is given by 


E 


o -A 0 2 


= 1- (l + A - " 2 ) 


-a 2 \- 1 / a 2 


A > 0. 


(29) 


The distribution of G(t) at time t > 0 is specified by 


E 


Gi(t) -X 2 G 2 (t) 


/ , n _q, 2 \- 1/«2 

(i + +a 2 J J ,0 


= 1 - 


< Si < 1,A 2 > 0, 
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and the transition probabilities for 0 < to < t\ < oo and (n\,y) € NqX R + are given by 


E 


sf 1 ^ ) e-^ G ^\G 1 (t 0 ) = n 1 ,G 2 (t 0 ) = y 
^- 1/ai (1 - Sl) i- Q1 ^ 


x exp 


= 11 - 
to 


1 - r + rf 1 - s iC“' 

to t o 


-l/ai' 


n i 


1-U + 


1 1 


-i 


ai 


1 — 1 /ai N 


+ A 2 \y 


In the sequel the symbol C^ will denote the law with initial condition p which may be a 
measure or a random variable, and we will write (•) for £(-|Z(n) 0). Moreover, we say 

that a process {0(f), f > 0} has a.s. constant paths with (random) value 0 if 

P(0(t) = 0 for all t € (0, oo)) = 1. 

We have the following global description of the reduced process under Condition ((HJ) and the 
assumption Z(0) = e2, where the symbol =>• means convergence in the respective space of 
cadlag functions endowed with Skorokhod topology. 


Theorem 4. Let Conditions ([2])-(|4]) and © be satisfied. 

(1) If m = [(1 — e _t )n] , t > 0, then 

jC^{Z(m, n), 0 <t < oo} =>• T( 0> i){(0, X(t)), 0 <t < oo}; 

(2) If m = n — th(n), t > 0, for h*(n) <C h(n) <C n, then 

C {n) | (zi(m,n), ^^^yyZ 2 (m,n)^ , 0 < t < ooj =>- £{(0,0 2 (*)),O < t < oo}, 

where the limiting process has a.s. constant paths with value 9 2 specified by V29\) : 

(3) If m = n — th*[n), t > 0, then 

C {n) | ^Zi(m,n), Q 2 ^ff w )) n )) > 0 < f < ooj £{ G (t),0 < * < °°}i 

For any Xi, X 2 > 0, any function h(n ) satisfying 1 -C h(n) <C h*(n) and m = n — th(n) 

£l "’ { n) ' ">) . 10 < * < °°} =*• £ ( 0 (‘). »<*<»}> 
where the limiting process has a.s. constant paths with value (9 1 , 0 2 ) specified by 


E 


g—Alfllg—A202 


= 1 - 



(w°‘ 



— 1/Q2 


Let d n and q n , n = 1,2,..., be positive functions such that d n q n = 1 and 

lim d n = oo and d n n~ £ = 0 for any e > 0 . 

71 —»• OO 

Then we have the following result under Condition Q: 

Theorem 5. Let Conditions and 0) be satisfied. Then 

( 1 ) 

{Z(tg* (n), n), 0 <t< oo} =>- £( 0i i){Y(t), 0 <t < oo}; 

(2) If ^ =02/(01(1 + 02)), then 

{Z(d n g* (n)n y , n), 0 < y < 1 -72} =>- £(<?i,o) {(^1(2/), 0), 0 <y < 1 - 72} 
where 9 1(-) has a.s. constant paths with value 9 \ specified by 


(30) 


E 


D 0i 


= 1 - (1 - ajVU+aa), 0 < S < 1. 


(31) 
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(3) If m = [((1 - e~ l ) + q n ) n\ , then 

£^ n ' ) {Z(m, n), 0 < t < oo } ==> £( 0 1 )O ){(Fi(i), 0 ), 0 < t < oo}, 

where the law of 9\ has been specified in ED- 

(4) If 1 <C h(n ) <C n and m = n — th(ii) then 

£ (n) { ( Qi{th{n)) Zl ^ ( ' m ’ n ' ) ) ’° < f < °°} t < °°}) 

where the limiting process has a.s. constant paths whose value p is specified by 


E 


0 -\p 


= l _ (1 _|_ / \-ai^-l/(«i(l+«2))_ 


For the case ot\ = a.i = 1, where 2g y 21 J (+i, S 2 ) = s 2 + s 1 , this result is established in 
Finally, we have the following result under Condition (11011 and the assumption Z(0) = e 2 : 

Theorem 6 . Let Conditions and m be satisfied. 

(1) If m = [(1 — e _t )n], then 

£("){Z(m, n), 0 < t < 00 } => £( 0 ,i){W (t), 0 < t < 00 }; 

(2) If n — m = th{n), where 1 <C h(n) <C n, then 


r(n) f ( Q2l{n) r, , , Q 2 {n) ^ ^ n „ , 

* 1 {n Ql (th(n)) Z ^ m ' n) ’ Q^m Z2{m ' n) ) ’° < * 


< 00 


£{Y(t), 0 < t < 00 }, 


where the limiting process has a.s. constant paths with value (v\, vf) specified by 


E 


e -Aivi e —A 2 i’2 


— 1 — ip (Ai, A 2 ). 


We end this section by the description of the distribution of the MRCA. More precisely, 
noticing that 

{Zi(m,n) + Z 2 {m,n) = 1} = {/3 n > m}, 

where the definition of j3 n has been given in ([Tj) , we deduce from Theorems [|]| 6 ] the following 
statements concerning the birth time j3 n and the type T n of the MRCA: 

Theorem 7. Let Conditions (Hlb|2P be satisfied. 

1) If Condition is fulfilled then for any a € ( 0 , 1 ) 

lim P (/3 n < an,T n = 2|Z(n) / 0) = o; 

n—>-oo 


2) If Condition is fulfilled, then for any t E (0,oo) and a E (0,1), 


lim Po 


< tg*(n),T n = 2 Z(n) / 0 ) = 


o 2 


( 1 -e 


-{l+a 2 )t/a2\ 


and 


1 + 0 2 

lim P 2 (s'*(n) < fin < an,T n = 1|Z(n) / 0) = 


3) If Condition lidl ) is fidfilled, then for any a E ( 0 , 1 ) 

1 09 

lim P 9 (/ 3 n < an|Z(n) / 0) = 1 — --(1 — a) — 

n->oo 1 + 012 


Moreover, the type of the MRCA satisfies: 

lim P 2 (T n = 2|Z(n) / 0) = 


1 + o 2 
«2 b a2 


1 + a 2 

_ a )((i+« 2 )&“ 2 -i)/a 2 _ 


(1 + 02 )6° 2 — 1 
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Remark 1. As a by-product of the proof of Theorem ® (see Equation m) we can obtain a 
simple expression for the death moment S n (2 ) of the last ancestor of type 2, 

5 n ( 2 ) := min{m < n : Z 2 (m,n ) = 0 } . 

Namely, under the conditions of Theorem® for any t E ( 0 ,oo) 

lim P 2 (5 n (2) < tg*(n)\Z(n) ± 0) = 1 - (1 + t)~ 1 / a \ 

n—>• oo 

Theorem [7] states that the law of the MRCA under Conditions @ and ( 1101 ) essentially differs 
from its law under Condition (JSJ), where it is almost surely of type 2 and with birth time 
uniformly distributed on the time interval [ 0 , n\. 


4. Auxiliary results 

In this section we list some known results and prove a number of auxiliary lemmas. The first 
statement, which will be useful at several occasions, is a direct consequence of the representation 
of regularly varying functions (see for example [7j, Ch.VIII, Section 9). 

Lemma 1. Let R(n ) be a function regularly varying at infinity with index —a <0. If n> m —> 
oo then R(m ) R(n) if and only if m <C n. 

We now reformulate Theorem 1 in (28], formula (11) in [ljj and Theorem 1 in (23] as a single 
statement. It will be an important tool in describing the local behavior for the reduced processes. 

Theorem 8 . Let Conditions (0)-^]) be satisfied. 

(1) If Condition (jSJ) holds, then 

Q2i(n) ~ Q2(n), (32) 

and for Ai > 0, A 2 > 0 


lim E [exp{—\iQi(h*(n))Zi(n) — \ 2 Q 2 (n)Z 2 (n)} |Z (n) / 0 ] 

n—>oo 


= 1- (l + (A}- Ql +A 2 )-“ 2 ) 

(2) If Condition Q holds, then 

Q 2 i{n) ~ Ql /{1+a2 \n)L 3 (n), Q 2 (n ) = o(Q 2 i(n)), 
where L^{n) is a slowly varying function as n —» oo, and 

Q^+“ 2 (n)L 2 (Q 21 (n))~ A 21 Q!(n), 

where L- 2 has been introduced in ([3]) . In addition, for all \\ > 0,A 2 > 0 

exp {—XiQi(n)Zi(n)} |Z(n) + ol = 1 - (1 + A^ 1 )" 1 /^ 1 ^ 2 )), 


-l /«2 


lim E 

n—too 


and 


1 


n-HKX> Q 2 (n) 


Eo 


1 — exp < —A 


Q 2 (n) 


n 


Z\{ri) - X 2 Q 2 (n)Z 2 (n) 


— 0(Ai, A 2 ). 


(3) If Condition (TTT71) holds, then 

Q 2 i(n) ~ bQ 2 (n), 

where b E (1, oo) has been defined in (fl9|) . In addition, for all Ai > 0, A 2 > 0 


n ->°o Q 2 i(n) 


E-; 


1 - exp <J Z I (n) - X 2 Q 21 (n)z 2 (n) 


n 


(33) 

(34) 

(35) 


(36) 

(37) 


if(Xi,X 2 ). (38) 
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To simplify notations we agree to write s k = for any vector s = (si,s 2 ) and integer¬ 

valued vector k = (k\,k 2 ). Besides, for vectors x = {x\,X2) and y = (2/1,2/2) we set x 0 y = 
(xiyi, X 2 IJ 2 ) an d let 1 = (1,1). Let us now present a way of expressing the probability generating 
function of the reduced process which will be repeatedly used in the proofs and will allow us to 
apply Theorem [7) To do this, we need to introduce some notation. Put, for i € {1, 2} 


Fi(n-,Si ) = E 


Z.i(n)\ 


Z(0) — e,; 


Qi{n\Si ) = 1 - Fi(n;si ), 


F 21 (n;s) = E[s z (fo|Z(0) = e 2 ], 
Q 2 i(u;s) = 1 - E 2 i(n;s). 


By conditioning on Z (m), we can write for m < n 


E 2 


g Z (m,n) 


F 21 (m\ 1 - (1 - si)Qi(n - m), 1 - (1 - s 2 )Q 21 (n - m)). 


Hence, setting Q 2 i(fc) = {Qi(k),Q 2 i(k )), k = 0,1,we get 


1 — E 2 


s z (m,n)|z( n ) ^ 0 


(Q 2 i(?r; 1 - (1 - s) 0 Q 2 i(n - m)))/Q 2 i(n). 


(39) 


Thus, to prove Theorems [1][3] it is necessary to find the limit of the right-hand side of (13911 under 
an appropriate choice of m, n and s. If n — m -A 00 it is equivalent to finding 


lim Q 2 i(m; 

m,n—>00 


e -(l-s)®Q 21 (n—i" 0 )/Q 21 („). 


(40) 


The end of this section deals with several expressions similar to (|39l) or (14011 . These results will 
be needed at several occasions in the remaining sections. 


Lemma 2. If m <C n, then, for A > 0 and i € {1, 2} 

Qi(m ; e ~ AQi(n) ) ~ A Qi(n) ~ Qi{m ; 1 - A Qi(n)). (41) 

Proof. Let h\(n) be the integer-valued function satisfying 

Qiihxin)) < 1 - e" AC Wfo < Qi(h\(n) - 1). 

We know by © that Qi{h\{n )) ~ Qi(h\(n ) — 1) and 

1 _ g-AQdfo „ AQj(n) ~ Q^nA" 1 ^), (42) 

implying h\{n) ~ nA _1 / Qi as n A 00 . Using the branching property we have 

Qj(m + h A (n)) = Qi(m,Fi(h\(n),0)) < Qi(m; e' AQi(n )) 

< Qi(m, Fi(h\(n ) - 1,0)) = Qi{m + h A (n) - 1). 
Since m <C h\(n), we get, again by © that, as n —> 00 

Qi (m;e~ AQi(n )^ ~ Qi{m + h A (n)) ~ Qi{h\(n )) ~ Qi^nA _1/ "^ ~ A Qi(n) 

proving the first equivalence in (1411) . The second equivalence follows from (1421) . ■ 


Lemma [3] deals with asymptotic results on the time scale h*[n ) specified in (Hill . 
Lemma 3. Let Conditions (Qj-Q) and © be satisfied and Ai,A 2 ,a be positive. Define 

s = 1 — \iQi(h*(n)). 

Then we have the following two equivalences for large n, 


n —1 


k=0 


A i-m 


A 


21 


Q 2 (n) 


ah*(n )—1 

X] Qi( fc ; s ) 

k=0 


\ 1 ^1 /■ \ 

(l-(l + aAr) 1 - 1 /ai )Q 2 (n), 


(43) 


(44) 




r\_/ 
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where we agree to understand ah*(n ) as [ah*{n)\. Moreover, 


lim 

n—> oo 


Q 2 \{ah*{n ); 1 


XiQi(h*(n)), 1 - A 2 Q 2 (n)) 
Qiiji) 


A } - " 1 (l - (1 + +A 2 . 


(45) 


Proof. We do not prove the first assertion in (1441) . as it may be checked in a similar way as 
the second one. As in the proof of the previous lemma, if h x {n) is an integer-valued function 
satisfying 

Qi(h\ 1 (n)) < 1 -s< Qi{h Xl (n) - 1 ), 

then 

Qi{k + h Xl (n)) < Qi{k ; s) < Qi(k + h Xl (n) - 1) (46) 

for every integer k. Equations (TiTfli and Q entail that, as n-> oo 


By (1461) we deduce the inequality 


Mn)~Ar ai fc*(n). 


(47) 


ah*(n)~ 1 ah*(n)+h\ 1 (n)-l 

0< ^2 ^2 Qi(k) ^ Qi{h Xl {i 

k =0 k=h\ 1 (n) 

Note that in view of (|5]) . as z —> oo 

oo 1 

~ Ql zQi(z) ~ -- z 1 ~ 1 ^ ai l\(z). 

f—■ ' 1 — cci 1 — ai 


- 1 ). 


k=z 


Hence we conclude that 


a h { n ) 1 oo oo 

E «i(*> = E Qi{k) — J2 

k = h \ 1 ( n ) k = h \ 1 ( n ) k = ah *( n )+ h \ 1 ( n ) 

~ 1 °2 (h\An)Qi(,h Xl (n)) - ( ah*{n ) + h Xl (n))Qi(ah*(n) + h Xl (n))^ 

~ ——^—A} _ai (l — (1 + oA“ 1 ) 1 ” 1/QlN )/i*(n)Q 1 (/i*(n)) 

1 — Oi V / 

~ ^A}- Ql (l - (1 + aA“ 1 ) 1 - 1/ai )Q 2 (n), 

where we have applied (|471) , ([5|) and (fill) . 

We now prove (j45j). By definition, for s = (si, s 2 ) E [0, l] 2 and fc € No, 

Q 2 i(k + 1; s) = Q 2 i(fc;s) - (52^ a2 (fe;s)L 2 (Q 2 i(fc;s) + (A 21 - p (k; s)) Qi(/c; si), 

where 

p(fc;s) = p(Fi(fc;si),F 2 i(A:;s)) —>■ 0 
as k —> 00 uniformly in s E [0,1] 2 - Thus, 

ah * (n) —1 

Q 2 i(ah*(n)] s) = Q 2 i( 0; s) - ^ 

k =0 

Let, for sufficiently large n 

si = 1 — XiQi(h*(n)) and s 2 = 1 — A 2 Q 2 (n). 

Using the inequality 1 — cic 2 < (1 — ci) + (1 — c 2 ) for ci,c 2 E [0,1], we get 

Q 2 i{k;s) < E 2 [Zi(k)\ (1 — s-J + 1 — s 2 =A 2 \k (1 - si) + 1 — s 2 . 


Qli a2 (k-s)L 2 (Q 2 i(k; s)) - (A 21 - p (fc; s)) Qi(fc; Sl ) 
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Adding (1321) and (fT2l) . we obtain for k < ah*(n ) and some C\ independent of n, 

Q2i(k] s) < XiA 2 ikQi(h*(n)) + A 2 Q 2 (n) < CiQ 2 (n). 

This, in view of the monotonicity of y 1 +a 2 L 2 (y) as y 0 and h*(n) <C n yields 

ah*(n)—l . . . . 

^ Qlt a 2 (t,s)L 2 (Q 2 i(fc;s)) < C 2 h*(n)Q 1 2 +a \n)L 2 (Q 2 (n)) ~ C 2 = o(Q 2 (n)), 


k=0 


a 2 n 


where C 2 is finite and independent of n and we have applied © to get the equivalence. Finally, 
recalling that 

ah* (n) — 1 

^ (A 21 - p (A:; s)) Q^k; Sl ) ~ A}“ Q1 (l - (1 + aA? 1 ) 1 - 17 * 1 ) Q 2 (n) 

k=0 

by (tBD and that 

Q21 (0; s) = 1 — s 2 = A 2 Q 2 (n) 
we get (1451) and end the proof of the lemma. ■ 


We now focus on time scales h(ri) larger than h*(n ), still under Condition ([5]). 


Lemma 4. Let Conditions and m be satisfied and h(n ) be an integer-valued function 

such that h*{n) <C h{n). Then, for any s E [0,1) 


lim sup E 2 

n ^°° m<n-h{n) 


1 - s Zl ^ m ’ n) \ Z (n) 0 


E 2 [1 - 

lim sup - L ——7—r-- = 0. 

n ->°° m<n-h(n) Q 2 I W 


Proof. Since Zi(m,n) is monotone increasing in m given n, it is sufficient to consider m = 
n — h{n). By conditioning on Z\ (m) we get 


Eo 




— Eo 


(Fi(n - m\ 0) + s (1 - F\(n - m; 0))) Zl(m) l > E 2 [pf l(m) (n - m; 0) 


Let C(k),k = 1,2,... be the total amount of type 1 particles at moment k produced by all 
particles of type 2 existing in the branching process at moment k — 1. Then, 


E 2 

Ef l(m) (n_m; 0) 

— e 2 

~m— 1 

F^ k \m — k; F\(n — m; 0)) 

— e 2 

~m— 1 

f[ F^ {k \n-k; 0) 




.k= 0 


. k=0 


By iteration, we also deduce for every k = 1,2,...: 


Eo [C(fc)] = E 2 [E 2 [C(k)\Z 2 (k - 1)]] = A 21 E 2 [Z 2 (k - 1)] = A 21 . 


Using the estimates above we obtain 


E 2 

l _ s z i (rn,n) 

< e 2 

1 

o' 

-if 

1 

S' 

_I 

1 

t-H 

1_ 

< e 2 

~m— 1 

Y C( k )Qi( n - k ) 




k =0 


. fc=0 


OO 

<A 21 Y Qi(k) = 0(h(n)Q 1 (h(n))), 

k=n—m 

where at the last stage we have used © and the fact that, as n -A 00 

Qi(n) = o (n~ 1 Q 2 (n)) = o (n _1 ~ e ) 

for some e > 0. Since h*(n) <C h(n) as n A 00, we have by Lemma [U (fTTT) and (1321) that 

h{n)Q\{h{n)) = o(h*(n)Qi(h*(n))) = o(Q 2 i(n )), 


proving the lemma. ■ 
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We now state asymptotic results on the time scale g*(n) specified by (US). 

Lemma 5. Let Conditions and (01 be satisfied and (si,s 2 ) be in [0, l] 2 . Then 

( 1 ) 

Q 2 {g*{n)) ~ a 2 A 2 iQi(n)g*(n); (48) 

(2) If m <C g*(n), then 


Eg 


1 — s' 


Zi(m,n) 

1 


(3) If g*(n) <C rre < re, then 


Eo 


1 — s 


Z 2 (m,n) 


= o(Q 2 i(n)); 


= o(Q 2 i(re)). 


(49) 


(50) 


Proof. First observe that according to (fTTD . < 7 * (re) <C n, and that (1451) is a direct consequence 
of (IS and (1351) . To prove the remaining statements note that by the branching property 


E 2 


1 — s 


Z\ (m,n) 

1 


< E 2 [Zi(m, n)] = E 2 [Zi(m)\ Q\ (re — m) = A 2 imQi (re - m). 


(51) 


Recalling the condition m<Cj* (re), Lemma [T] and (1481) , we conclude that 

1^1 ( i *( i/D i */ \\ */ \ Q 2 {g (t )) Q2l(n) 

rreQi (re - m) < g (n)Q 1 (re - g (re)) ~ g (n)Q 1 (re) ~~- - —. 

ol 2 A 2 1 a 2 A 2 1 

This proves (US- To justify (15Uj) it is sufficient to observe that, for g* (re) <C m, 


Eg 


1 — s 


Z2(m,n) 


< Q 2 (m) < Q 2 (g*(n)) ~ Q 2 i(n). 


This ends the proof of the lemma. 


To end this section we mention a classical result on branching processes. Recall Definition [T] 
and define {M = (M 1 (t),M 2 (t)),t > 0} = {re, (gi(si),m), (g 2 i (s),/x 2 )}. Let 

/i(t;si) = E[sf lW |M(0) = ei ], / 2 i(f;s) = e[s m «|M( 0) = e 2 . 

Then we have the following statement (see, for instance, [12], Ch. IV, Section 3 or [T] , Ch. V, 
Section 7). 

Lemma 6 . The pair of functions (/i,/ 2 i) is the unique solution of 

{ df2 dt ’ S) = ^2 (521 (/1 {t m , si) , /21 (i; s)) — / 2 i (f; s)), 

\ 9fl m Sl) = vi{gi{fi{t;si)) - /i(t;«i)) 
with initial conditions f\ ( 0 ;si) = si and / 2 i( 0 ;s) = s 2 . 

5. Proof of Theorem |T| 

Throughout this section we assume Conditions ©-(HI) and (0), and that (si,s 2 ) belongs to 

( 0 , 1 ) 2 . 

Point (0): Let m <C re. In view of (0) and Lemma [4] 

Z 2 (m,n) | 


0 < E 
Thus, 


|Z(re) 0 


- E 


s z ( m , n )| Z ( re ) ^0 


E 


s z ( m , n )| Z ( re ) ^0 


= 1 - E 


1 — s 


< E 


Z 2 (m,n) 


1 — sf 1 ( m ’ n )|Z(re) 7 ^ 0 =o(l). 
/Q 2 i(n) + o(l). 
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Further, using the branching property and Lemma [2] we get 


E 


1 — s, 


Z2(m,n) 


= Q 2 (m; 1 — (1 — s 2 ) Q 2 {n — m)) ~ (1 - s 2 ) Q 2 (n — m) ~ (1 — s 2 ) Q 2 (n). 


This, on account of (1321) implies (fl3l) . 

Point (1): Let a be in (0,1). As in the previous case, by 0 and LemmaQ] as n —>• oo, 


E 


z (an,n) | 2 ( n ) ^ 0 1 = E \ s ^ an ’ n) \Z (n) / oj + o(l). 


Moreover, in view of J7D and (152|) . as n —> oo 


1 - E 


Z^ian^n) 


Z(n) + 0 


Q -2 (an -1 - (1 - s 2 )Q 2 i((1 ~ a)n))) 

Q2i(n ) 

Q 2 (an, e -{l-S 2 )W{l-a)) 1/a *Q 2 (.an)) 


a 1 /°‘ 2 Q 2 (an) 

= (a + (1 — a)(l — s 2 ) -Q2 ) -1/a2 


l/«2\ —«2\ —1/CK2 


~ a 


This completes the proof of point (1). 

Point (2): Let n — rn = h{n), where h(n) is an integer-valued function such that h*(n ) <C 
h{n) <C n. Similarly to the previous point we have by ([HD and Lemma 0] 


0 < E 


Z2(m,n) | 


Z(ra) / ol — E [s z(m ’ n) |Z(n) / ol = o(l) 


as n —» oo. Now following the line of proving point (1) and letting 


1 — s 2 = 1 — exp < —A 2 


Q 2 (n) 


A 2 - 


Q 2 (n) 


Q 2 i(h(n)) J Q 2 i(h(n )) 


we get, using again (ED, 


lim E 

n—>oo 


l — g Z(n—/i(n),ra) 


Z(n) + 0 


= lim 

n—>• oo 

= lim 

n—>• oo 


Q 2 (n - h(n); 1 - (1 - s 2 )Q 2 i(/i(n))) 

Q2i(n) 

Q 2 (n — /r(n); 1 — A 2 Q 2 {n)) 


Q2(jl) 


= (I + V”) 




Point (3): Let t be positive and m = n — th*{n). We need to find the limit of the right hand 
side of (1391) with 


1 — s 2 = exp < —A 2 


Q2(n) 


A 2 - 


Q 2 {n) 


Q 2 \{th*{n)) J Q 2 i(th*(n )) 


From (1331) we know that 


Hm Q 21 (m;e- x ^ hf ^-e- x ^^) = / + 


-l/«2 

_ : : 12 = 1 1 m t X 1 -"! rid 

Q 2 i(m.) 

Hence, by taking a fixed si E [0,1) and adding (0 we deduce for m = n — th*(n ) that 
Q 2 \(n — th*(n ); 1 — (1 — s) <g> Q 2 i (th*(n))) 


lim 

n—>■ oo 


= lim 

n—>• oo 


Q 2 i(ra) 

Q 2 i(n — th*(n ); 1 — t _1 / ai (l — si)< 5 i(h*(n)), 1 — A 2 Q 2 (n)) 

Q2i(n) 

= (1 + ( t 1 " 1 /" 1 (1 - si ) 1 ” 01 + A 2 )- Q2 )- 1/a2 , 


which is the statement of point (3). 
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Point ( 4 ): Take n — m = h(ri) with 1 <C h(n ) <C h*(n), 


, > Qi(A») 


and s 2 = exp < — A 2 


Qiin) 


Using (fTTTTI) we get 

Q 2 i(n- h(n); 1 - (1 - s) ® Q 21 (/i(re))) 

lim --——- = lim 


Q2i{h{n)) 
Q2i(n — h(n ); 


Q2i{n) 


Q 2 \(n - h(n)) 
= (l + (A}“ ai +A 2 )- Q2 ) _1/ “ 2 . 


This ends the proof of Theorem |T] 

6. Proof of Theorem \2\ 

Throughout this section we assume that Conditions ©-gD and (J21) are in force and that 
(si, s 2 ) belongs to ( 0 , l) 2 . 

Point ( 0 ): Let m -C < 7 * (re) with g*(n) specified by (USD. Applying g 9 |) we get, as re —» 00 


0 < E 


1 - s 


Z(m,n) 


- E 


1 — s 


Z 2 (m,n) 


< E 


1 — S 


Z\{m,n) 


= o(Q 2 i(n)). 


Further, recalling (11611 and Lemma [2] we see that 


E 


1 — s 


Z 2 {m,n) 


= Q 2 (m; 1 - (1 - s 2 )Q 2 i(re - m))) 

~ Q 2 (m; 1 - (1 - s 2 )Q 2 (g*(n))) ~ (1 - s 2 )Q 2 (g*(n )) ~ (1 - s 2 )Q 2 i(re). 


Thus, 


_ , % ttVi „Z(m,n)l rpri Z 2 (m,n)i 

T XT'Tl Zl(jn^ri) I ry r \ / ~-i y E l S J E[1 S 2 J 

lim E l - s 2 ZIn) ^ 0J = lim -——- = lim - — " -= 1 - s 2 , 

n—yoo n—yoo Q 2 i{n) n^roo Q 2 l{n) 

which proves point ( 0 ). 

Point (1): In view of Q and (1161) . for any t > 0 

Q 2 i(n - tg*(n )) ~ Q 2 i(n) ~ Q 2 (g*{n)) ~ t 1/a2 Q 2 (tg* (re)) 
and by (l48ll and © 

n * *,„x\ r\ Q 2 {g*{n)) t 1 / 012+1 Q 2 {tg*(n)) 

Qi(n - tg re ~ Qi(n) - - -- ~- - - . 

a 2 A 21 5 (re) a 2 A 21 tg*(n ) 

Then, by using m , csd and (I36|) with re replaced by tg*(n) we obtain 

E 2 m _ s z(t 9 *(n),n)i Z (re) / 0 ] = 1 ~ (1 ~ s) (8 Q 2 i(n - tg*(n))) 

2 Q 2 i(n) 


E 2 


1 — exp j — (1 — si 


p/ Q2 + i Q 2 (tg*(n)) 
a 2 A 2 i tg*(n) 


l (tg* (re)) - (1 - S2)t 1 / a2 Q2(tg*(n))Z 2 (tg*(n))} 


t 1 /°‘ 2 Q 2 (tg*(n )) 


*- 1/aa 0((l - (1 - S2)t 1/a *) 


a 2 A 2 i 


justifying the statement of point ( 1 ). 

Point (3): If g*{n ) <m<n then, according to (15011 


E-; 


m,n) _ Q Z 2 (m,n)\ ^ N _ 0 ^ 2 (m,n)l = 0 (Q 21 (re)) 


Hence we obtain the same limit for (11811 independently of the choice of s 2 . We choose s 2 satisfying 
(1 - s 2 )Q 21 ((1 - a)re) ~ Q 21 ^^ ~ (1 - si) 1/(1+ " 2) (l - a) _ 1 /( “ l( 1 + “ 2 )) Q 2 i(re), 
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where we have applied (1M1) . Then, using the branching property, © and (ITU) , we get 


Q 21 (an; 1 - (1 - s) <g) Q 2i ((l - a)n)) Q 2 i(an; 1 - Q 2 i((1 - a)n/(l - s i)" 1 )) 


Q2i(n) Q2i(n 

_ Q 2 i(an + (1 - a)n/(l - si) Ql )) 

Q2i(n) 


a + 


(1 — a) \ -l/(ai(l+« 2 )) 


( 1 -Sl )" 1 


Adding ([39]) . this proves point (3). 

Point (2): We know from point (1) that for every t > 0, 

lim P (Z 2 (tg*(n),n) = 0|Z(n) ^ 0) = 1 - t~ 1/a2 ^{0, t 1/a2 ). 

n—>oo 

Setting q(x) = (j>( 0,x)/x, x > 0, we deduce by applying (fT5l) that q(x) satisfies 

q'(x) = —x ot2 ~ 1 q 1+a2 (x), 

with initial condition q( 0) = lim^o q( x ) = 1- This equation has an explicit solution and we 
obtain 

lim P 2 (Z 2 (tg*{n),n) = 0 |Z (n) ± 0) = 1 - q{t 1/a2 ) = 1 - (1 + t)~ l/a2 . (52) 

n—>• oo 

In particular, 
and for rn g* (n), 


lim lim P 2 (-Z 2 (t<?*(n), n) = 0 |Z(n) / 0 ) = 1 , 

t—>oo n —>-oo 


lim P 2 (Z 2 (m, n) = 0|Z(n) / 0) = 1. 

m,n—^ oo 


(53) 

(54) 


We now take r(x) = 4>(x,0)/x 1 ^ 1+a2 \x > 0. Then from point (1) we have 


lim E 2 


Z 1 ((a 2 A 21 x/(l-s 1 ))°‘2/(l+<*2)g*( n ),n) 


Z(n) / 0 


= 1 - 


/ 1 — S]\ l/(«2 + l) 

Va 2 A 2 i ) 


r(x). 


As Z\(m, n) is non-decreasing with respect to m and we are looking for non-negative quantities, 
it follows that r is a non-negative, non-decreasing and bounded function. Hence it admits a 
non-negative limit at infinity. Moreover from (1531) we see that 


lim lim E 

t —>oo n—too 


QZl((a 2 A2id“ 2/(1+Q2) 9*( n ). rl )|2;(n) / 0 


= 0 


in view of Zi(m,n) + Z 2 (m,n) > 1. Hence lim x ^. 00 r(x) = (a 2 A 2 i) 1 ^ Q2+1 ' ) and 


lim lim Efs 

t —>oo n—>oo 


Zi(tg*(n),n) 


1 


Z(n) / 0] = 1-(1-Si) 1/(1+Q2) . 


From point (3) we also know that 

lim lim E[sf i( “"’ , ^|Z(n) / 0] = 1 — (1 — si) 1 ^ 1-1 " 02 ^. 

a—>0 n —>oo 

As Z\ (rn, n) is non-decreasing with respect to m, we end the proof of point (2) with (1541) . 
Point (4): The proof is similar to that of point (3) and we give less details. In particular, 
the choice of s 2 has no influence on the limit, as there is no type 2 individuals anymore in the 
reduced process. Let n — m = h(n) with 1 <C h(n) <C n. Set 

i _ i „ f ^ Qi( n ) 1 x Qi ( n ) 

51 XP { Al g!(Mn)))~ Xl Qi(Hn)) 

and choose s 2 such that 

(l- 8 2 )Q 2 i{h{n))~Q 21 (\l a 'n). 













REDUCED TWO-TYPE BRANCHING PROCESSES WITH INFINITE VARIANCE 


17 


Then by (1391) . the branching property, © and (f3~iT) we conclude that 

Q2i(n - h(n ); 1 — (1 - s) <g> Q 21 (h(n))) 


E[1-s z M| Z („)^o] = 


Q 2 i(n) 


Q 2 i(n - h(n )-1 - Q 2 i(A- Ql n)) _ Q 2 i ((1 + A ~ ai )n - h(n)) 


Q-2i(n) 

which ends the proof of Theorem [2j 


Q2i{n) 

7. Proof of Theorem \ 3 \ 


~(1 + A 1 “ 1 ) “l(l+“2). 


Throughout this section we assume that Conditions (12j) ~ 

(si,s 2 ) belongs to (0, l) 2 . 

Point ( 0 ): Using m, Condition (11011 and (1371) we obtain for m <C n 


and CED are satisfied and that 


E 


1 — s 


1 


< A-2\mQi (n - m) <C nQi (n) ~ crQ 2 {n) ~ ab Q 2 i(n), 


implying 


0 < E 


l — s Z(m,n) 


- E 


1 — s 


Z2(m,n) 


< E 


1 — S 


Z\{m,n) 


= o(Q-2i(n)). 


Further, by Lemma [2] and (1371) 


E 


1 — s 


Z2{m,n) 


= Q2 ( m ; 1 - (1 - s 2 )Q 2 i(n - m )) ~ (1 - s 2 )Q 2 i(n), 


giving 

Point ( 1 ): The result is a direct consequence of (|38|) . (1391) . and of the equivalence Q\{n) ~ 
&Q 2 i(n)/b under Condition (flOl) . 

Point ( 2 ): Again it is a direct consequence of expressions (1381) and (1391) . As it is very similar 
to the proofs of points (3) and (4) of Theorem [2j we do not give the details. 

8. Convergence of finite dimensional distributions 

In this section we study the limiting behavior of the finite-dimensional distributions of the 
properly scaled reduced process {Z(m, n),0 <m< n}. To simplify notation we let 


s) = E s z(m ’ n) |Z(n) / 0, Z(0) = e, : 


i = 1.2 


and, given 0 < k$ < k\ < ... < k p < n set k = (/cq, k ±,..., k p ). Putting S 1 = (su, s 2 i ) denote 


f(k 0 ,k i,...,fc p ,ri), 


TV' v u,"'i ) ...) n 'P!'0 /o Cl _ hk.n) / o Cl _ TP 

J 2 ,...,Opj — J 2 foi, ..., Op) — Cj 


ns; 


Z(fci.n) 


J 


{k °*' = ii k '” ) (sii,..,si p ) = E 


LZ=l 
P 


n 

L=i 


Z\(ki,n) 


'll 


Z(k 0 ,n) = e 2 
Z(/c 0 ,n) = ei 


and 


j( k » n )(S 1; . , Sp) = (j^ n \su,..,3i P ), jf’ n) (S!,...,S p )) . 

The next statement is a simple observation following from Corollary 2 in 
Lemma 7 . For any 0 < ko < k± < ... < k p < n we have 

jf’ n) (Si,...,S p ) = jf 0 ’ fe l’ n) (Si <g> j( fcl ’ fc2 --^’ n )( S2 ,..., Sp )) 

= jfi-fco,n-fco) ^ Sl g, j(fc 2 -fcl,n-fci) ^g 2 g, f ^S p _i <g> j(*p-*P-i.»-*P-i)(S p )) ... 
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In particular, if k = ( 0 , k\, k 2 ), then 

{S l ,S 2 ) = J 2 (fcl,n) (Si® J (fc2 " fcl ’ n - fcl) (S 2 )). 

Remark 2 . It follows from Lemma [?] that the process {Z(to, n), 0 < to < n| Z(n) 7 ^ 0 } is, /or a 
given n, an inhomogeneous discrete-time Markov branching process in m. 


Observe that for k = (&o, k\) 


E 


i Z(kl,n) \Z(k 0 ,n) = (r 1 , r 2 ) 


= •/ 


r (k,n) 


ri 


Jr 


(k,n) 


(si,s 2 


r 2 


Hence to investigate this quantity we can use the equalities 


=1 - 


Qi(&i - fc 0 ,1 - (1 - si)Qi(n - fci)) 


and 


J. 




(s) = 1 - 


Qi{n - hi) 

Q 2 i(ki - k 0 ] 1 - (1 - s) <g> Q 2 i(n - ki)) 


(55) 

(56) 

(57) 


2 " Q 2 i(n-k 0 ) 

We are now able to prove convergence of the finite-dimensional distributions of the prelimit¬ 
ing processes to the corresponding finite-dimensional distributions of the limiting processes in 
Theorems Q] to EJ 

8 . 1 . Convergence of finite-dimensional distributions in Theorem [ 4 ]. 

Point ( 1 ): For t > 0 introduce the function 


f 2 i(t;s 2 ) = 1 - q - e 1 + e 4 (1 — s 2 ) 02 
whose first order derivative with respect to t is 
! = -L(e-<_e -* {1 
,(*) a 


= lim Eo 


Z((l-e 4 )n,n)| Z ^ ^ Q 
- 1 - 1 /ai 


1 -e " 4 + e~ t (l -s 2 )-" 2 ) 


dt a 2 

Recalling the definition of gf^ ] in ()23l) we get 

1 

OL 2 

' o-Rl - _ (1 _ e ~t + e ~t^ _ S2 ^-a 2 )-l/a 2 ' 

df 2 i(t;s 2 ) 


gi X \f 2 i(t;s 2 )) - f 2 i(t; s 2 ) = ^-((1 - f 2 i(t; s 2 )) 1+a2 - 1 + (1 + a 2 )/ 2 i(i; a 2 )) - f 2 i{t;s 2 ) 


= — ((1 — e _t + e _i (l — s 2 ) _ “ 2 )' 
a 2 


-"21 _ 


dt 


= —(1 - e~* + e-*(l - S 2 )-"2)-i-i/"2 _ (1 _ e -« + e -*(i - S2 ) 

a 2 V 

Applying Lemma [B] we conclude that for s\, s 2 E [0,1] 

lim j(°.( 1 -e“ t) ’M0( s ) = / 2i (t; s 2 ) = E L YW |X(0) = 1 

n —/oo L 

where A'(-) has been specified in Definition [TJ This, along with Theorem [1](1) proves convergence 
of one-dimensional distributions of the prelimiting process to that of X(-). 

In view of Lemma [Tj to check the needed convergence of arbitrary finite-dimensional distri¬ 
butions it is sufficient to consider the case p = 2. Let 

k\ = [(1 — e _tl )n] < [(1 — e~ t2 )n] = k 2 , 0 < t\ < t 2 . 

By the previous estimates, Lemma [71 Theorem [U(l) and the uniform continuity in t\ and t 2 of 
the functions involved, we get 

lim J 2 (k ’ n) (Si,S 2 ) = lim J 2 (Al>ri) (Si<g> j( A2 ’ n - fcl) (S 2 )) = f 21 (ty,s 21 lim J^ n ~ kl \ S 2 )) 

= /21 (tr, S21/21 (t2 — tr, s 22 )) = E = 1 
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as desired. 

Point ( 2 ): Convergence of one-dimensional distributions has been established in Theorem[l|2). 
In particular, there is no type 1 individuals in the limit process. As before, to check the needed 
convergence of arbitrary finite-dimensional distributions it is sufficient to consider the case p = 2 
and to calculate the transition density of the limiting process for type 2 individuals. Let h be 
an integer valued function such that h*(n) <C h(n ) <C n, and ko = n — toh(n) < k\ = n — t\h(n). 
Take s 2 satisfying, for a positive A 


1 — S 2 = 1 — exp 


x Q 2 (n ) } ^ A Qa in) 

Q 2 i(t 2 h(n)) j Q 2 i(t 2 h(n))' 


Let, further, 


Z 2 (k 1 ,n) = y 


Q 2 i{tih(n)) 

Q2{n ) 


(recall that we agree to consider the equalities of the form Z (x, n) = y as Z ([#] , n) = [y]). Then 
from (1561) and the fact that there is no type 1 individuals we get for sj in [0,1], 


i - i 2 (k ’” ) (s) 


1 — Eo 


o-^Q2(n)Z2({to-ti)h(n)) 


At? 2 (n) 


Q 2 i(t 0 h(n)) 

where we applied Lemma [2j Hence we deduce that 


Q 2 i{toh(n)) ’ 


log E 


Q2(n) 


Q2i(t 0 h(n)) ? (k,n)/ x x 

= y — rT7 ~\— logJ 2 ( s ) ~ ~ X V- 

Q2(n) 


Point ( 3 ): Once again, it is sufficient to calculate the transition density of the limit process 
only. Let ko = n — toh*(n) < k\ = n — t\h*(n ), where h*(n) has been defined in (fT2l) . Applying 
first (1561) we get 


i- jj k ’”W 


1 — Ei 


D -(l-s 1 )Qi(t 1 h*(n))Zi((t 0 -t 1 )h*(n)) 


1 — Ei 


Qi{t 0 h*(n)) 


({to - H)Ao) 1/ai Qi((io - ti)h*(n)) 


t() \ aj 


to — t\ 


1 + 


t\ 


to — t\ 




tl , 1 1 


to to 


Introduce now 


s 2 = exp < —A 2 


Q2(n) 


‘Q 2 i(tih*(n)) J ' 


Then, by Markov property and ()39l) we get 

l-jf’ n) (s) = 1 -E 2 (t 0 h*(n)) + 0 

Chi((to - ti) h*(n)\ 1 - (1 - s) ® Q 2 i(tih*(n))) 

Q 2 i(toh*(n)) 

Q 2 i((to - tl) h*(n ); 1 - (1 - si)t^ 1/Ql Qi(/i*(n)), 1 - A 2 Q 2 (n)) 


Q 2 i(toh*(n)) 


Q 2 (n) 


Q 2 i(toh*(n)) 


l—n 




to 


1 - 2 -, 


+ A 2 
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where for the last equivalence we applied (|45l) with a = to — t\ and Ai = (1 — si)^ 1 ^ ai . Hence, 


E 


jf (fe ' n) e- A2 02 i?*i ( fc*(«)) Z2{kl ’ n) \Z(k 0 , n) = y ^ x ° h *^ 

Q2{n) 


e 2 


exp 


- (a- s^h 1 - 1 ^ 1 (i - S1 ry 1/ai ^j + a 2 ) 


y 


This proves the desired convergence of the finite-dimensional distribution of the prelimiting 
process to those of G(-). 

Point ( 4 ): Let h be an integer valued function such that 1 -C h(n ) <C h*(n), and Ai,A 2 be 
positive numbers. Selecting ko = n — toh(n) < n — t\h(n) = k±, we get by using (1571) 


1 _ j( k 0,kl,n) A 1 Ql( tl h(n)) ) e A 2 Q 2 l(Hh(n))^ 


Q 2( n ) 


1 — E? 


exp {-A 2 Q 2 (n)Z 2 ((t 0 - h)h(n))} 


E 


21 


+ 


Q 2 i(t 0 h(n)) 

e -\ 2 Q 2 (n)Z 2 ((t 0 -ti)h(n)) f y _ e -AiQi(/i*(n))Zi((t 0 -ti)/i(n)) 


Q 2 l(t 0 h(n)) 

By applying twice Lemma [2] as h(n) <C h*(n) and h(n) <C n we get that the first summand is 
equivalent to 

A2<32(^) 

Q 2 i(toh{n))' 

and the second summand is smaller than a term equivalent to 


AiQi(fr») 

Q 2 i(t 0 h(n)) 


= o 


Q2(n) 


This yields 

E 2 


. Zi(ki,n)Qi{h*(n)) Z 2 (ki,n)Q 2 (n 
exp — Ai- ^ ^ , / __ \ \ -A 2 


Qi(tih{ri)) 

In the same way we deduce that 

Z 1 (ki,n)Qi(h*(n)) 


Q2i{hh(n )) 


Q 2 i(t 0 h(n)) 

Z(/c 0 ,n) = y 2 


Q 2 i(t 0 h(n)) 


Q2{n) 


1 — E 5 


exp < — Ai- 


Qi(hh(n)) 


\ Qi{t 0 h(n)) 
m 0 ,n) = yi Qi(ft . (n)) e! 


-e 2 


= — A u/i 


0 —\2V2 


The two last equivalences justify the claimed convergence of finite-dimensional distributions. 


8.2. Convergence of finite-dimensional distributions in Theorem O 
Point ( 1 ): Applying Theorem [2]^ 1 ) we can consider the function 

A +1 

ol 2 A 2 1 


/ 21 (f;s) := lim E 2 [ s z(tff* ("),*») i Z ( n ) / ol = 1 - t 

n —>-00 L 

whose first order derivative with respect to t is 


1 - Si) 


(1 - s 2 )t“a), 


df 2 i(t;s) 

dt 


1 .-±-1 

= - 1 “2 

OL 2 

1 1 


Ai, A 2 ) — (1 + «2 )Ai^y-(Ai, A 2 ) — A 2 —^-(Ai, A 2 

UA\ OA2 


OL2 


(t “2 y +a 2 (Ai,A 2 )-i + Sl ), 
where we have applied (fT5l) and used the notations 

Ai = (1 - si)t 1/,Q2+1 /(a 2 ^ 2 i) and A 2 = (1 — s 2 )t°v . 
Now taking f\ (i; si) = si and recalling the definition of \ ^ in (1241) we obtain 


g { 21 \sij 2 i(t]s)) -/ 21 (t;s) = 


«2 df 2 i(t;s) 
1 + a 2 dt 
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Applying Lemma [ 6 ] yields that /21 is the generating function of the process Y specified in 
Definition |TJ Hence the desired convergence of finite-dimensional distributions follows by the 
same arguments as before. 

Point (2): We take g*(n ) <C k^ < k\ <C n. From Theorem [2^2) we know that there is no more 
type 2 individuals in the limit process. Applying (|56l) and Lemma [2] we get 

i(k ,n), x , Qi(ki - k 0 , 1 - (1 - si)Qi(n - ki)) (1 - si)Qi(n - ki) 

J 1 (si) = l-777-77-~1-777-77-~»i- 

Qi{n-k 0 ) Qi(n-k 0 ) 

Point (3): Again as there is no type 2 individuals it is sufficient to study the dynamics of type 
1 individuals. Applying Theorem [2]^ 3) we can consider the function 

fi(t; s\) = lim E 2 [s z (( 1 - e_ ‘) n ’ n )|Z(n) ^ ol = 1 - ((1 - e" 4 ) + e“ t (l - si)-^ 1 )- 1/ “ l(1+Q!2 ) , 

n—too L J 

whose first order derivative with respect to t is 
dfi(t;si) 1 


dt «i(l + 0 : 2 ) 

Recalling the definition of * in (1251) we obtain 

1 


(e-* - e _t (l - Sl )- Q1 ) (1 - e~ f + e~\l - si)- ai )~ 1_1/(ai(1+a2)) 


9i } (/i(i;si)) - a x + Q!2 ^ ((!■ - /i(t;si)) 1+Ql(1+Q2) - 1 + /i(t;«i)) = 


dfi(t-,si) 
dt ' 


Applying Lemma E] yields that f \ is the generating function of the process V{-) specified in 
Definition QJ3). Hence the desired convergence of finite-dimensional distributions follows. 
Point ( 4 ): Let h be an integer-valued function satisfying 1 <C h(n) <C n, and ko = n — toh(n ) < 
k\ = n — t\h(n). The one dimensional distribution of the limit process is given by Theorem[2{4). 
Applying Lemma [2] and (1561) yields 


1 —E 
Hence we deduce 

as required. 


Xl Q®tihU Zl ( kl ’ n) \Z(k 0 ,n) = ei 


1 — Ei 


e —^iQi(n)Zi((to—ti)h(n)) 


Qi(t 0 h(n)) 


Qi(n) 

1 Qi(t 0 h(n))' 


-Ai t, ? l( u) v; Zi(ki ,, 
e Qi(*i h («)) v J \Z(ko, n) = y 


Qi(x 0 h(n)) 

Qi(n) 


ei 


0 -^iy 


8.3. Convergence of finite-dimensional distributions in Theorem [6l 
Point ( 1 ): Consider the function 


1 - /i(t; si) = 1 - lim Ei 
n—>00 


T 


Z(n) + 0 


1 — Ei 


= lim ■ 

h\ —^00 


„-(l-si)(e t -l) 1 /“lQ 1 (fc 1 )Z 1 (fc 1 ) 


— (l — e + e (1 — si) 


— Ot\ \ a l 


(1 — e-*)i / ai Qi( kl ) 

where we applied (1561) with ko = 0 and k\ = [(1 — e _t )n], and ©• Then following the proof of 
Theorem 2)^1) we obtain 


g[ U \h(t-,si)) - fi(t; si) = 


dt ' 


Now applying Theorem [3^1) we can consider the function 

/ 21 (f;s) := lim E 2 |>((i-e->,n)| Z(n) + ol = 1 - (1 - A 2 ), 

n—> 00 L 

where, for the sake of readability we have introduced 

Ai = (1 - s 1 )(e t - l) 1/ai and A 2 = (1 - s 2 )(e* - l) 1/c * 2 . 
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Since 1 + a 2 1 = a 1 1 under Assumption m > we get 

= (1 - e" 4 ) 


9/21 (t;s) n _ p -U-l/a 2 ~l(^_ 


dt 

(1 - e -t)-l/aa-l 
«2 


t 1 \ \ d i’ X \ A 2 9^ T 

^(Ai, A 2 ) A 2 ) tt- (Ai, A 2 


a 2 


or 0 Ai 


a 2 <9A 2 


e - 4 - 1)^(A!,A 2 ) + 6 " 2 t// + “ 2 (Ai, A 2 ) - a 2 A 12 Ai (1 + ( -Ai 


where we have applied (12011 . Further, recalling the definitions of and k in 
see that 

3 2 i ) (/i(Mi),/ 2 i(£;s)) - / 2 i(i;s) 

1 /ol4i 2 „ , . 6" 2 


ai^ — l/oi^ 

and 


h 


we 


Observing that 

1 b a2 (l + a 2 ) 

K 


fi{t\si) + — ((1 - / 2 i(t;s )) 1+ “ 2 - 1 + (1 + a 2 )/ 21 (t;s))) - f 2 i(t- s). 

6" 2 (1 + a 2 ) 


- 1 = 


«2 


a 2 (1 + a 2 )b a2 — 1 a 2 

and recalling the definition of b in (11911 leads to 
,(W)/ 


- 1 = 


1 


1 


(1 + a 2 )b a2 — 1 a 2 /t 


021 si), / 2 i(t; s)) - / 21 (i;s) 

^(> 2 (i - / 2 i(i;s )) 1+ “ 2 - - / l( t; Sl )) - 1 + / 21 (i;s)) 

cra 2 Ai 2 


a 2 n 

(1 - e-*)" 1 " 1 /^ 


o 2 k 


1+02 


(Ai, A 2 )— 




l+i-“Z7(l —s l) 


where we have used again the equality 1 + a 2 1 = aq i . This, finally, yields 


1 — e t 
-l 


—ai \ a l 


— (1—e )V’(Ai,A 2 ) , 


9 2 i\fi{t-,si)J 2 i{t-,s)) - / 2 i(t; s) = - 


1 9/ 2 i(t; s) 


k <9t 

Applying Lemma [ 6 ] we conclude that / 2 i is the generating function of the process W (•) specified 
in Definition [DA). 

Point (2): Let 1 <C h(n ) -C n and ko = n — toh(n) < k\ = n — ti/i(n) for 0 < U < to- Choose 
Ai,A 2 > 0. Using (f56l) . the fact that h(n ) 1 and recalling Lemma [2] we see that 

° l(n) N Qi((/o - ti)h{n), 1 - AiQi(n)) AiQi(n) 


1 _ j( k 0M,n) A l Ql (/ i r(n))^ 


Qi(x 0 /i(?r)) 


Hence 


logE 

Then from ()57H we obtain 

V Qi (") 


e Al QidiM")) Zl(fcl,n) |Z(fc 0 , n) = y 


Qi{t 0 h(n)) 

Qi(n) 


ei 


Qi(x 0 h(n ))' 
-Ai y. 


. /, , \ / Ai y-i ini AoQo(ti) \ 

1 _j{ k 0M,n) f £ ~ Q l(tl h(n)) ^ e ~ Q 2 (t 1 h(u)) \ 


1 — Ed 


?—XiQi(n)Zi((to—ti)h(n)) p —\ 2 Q 2 i(n)Z 2 ((to—ti)h(n)) 


1 — Ed 


o—X2Q2i(n)Z2((to—ti)h(n)) 


Ed 


+ - 


Q 2 i(t 0 h(n)) 

e —*. 2 Q 2 i(n)Z 2 ((to—ti)h(n))f ^ _ g—AiQi (n)Z\((to—t\)h(n)) 


Q 2 i(t 0 h(n)) Q 2 i(t 0 h(n)) 

To control the first term we apply Lemma [2] and Equivalence (13711 : 


1 — Ed 


o~^ 2Q21 {n)Z 2 ((to-t 1 )h(n)) 


Q 2 i(t 0 h(n)) 


A 2 Q 2 i(n) 

Q 2 i(t 0 h(n)) ’ 
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The second term is non-negative and bounded by 


E 2 

l _ g-AiQi(n)Zi((t 0 -ii)/i(n)) 

Ai(to - ti)h(n)Qi(n) _ 

( A 2 Q 21 {'a') \ 

Q 2 i(t, 0 h(n)) Q 2 i(t 0 h(n)) 

^Q2i(t 0 h(n))J 


where the equality is a direct consequence of Condition ([TUI) and Equivalence (f37l) . This ends 
the proof of point ( 2 ). 


9. Tightness 

In this section we complete the proof of Theorems SK] by checking tightness. A key tool will be 
a slightly simplified version of Theorem 6.5.4 in |9j giving a convergence criterion in Skorokhod 
topology for a class of Markov processes which we formulate to simplify references. Let ai < a 2 
be in R + and 

{K n (u) = (K ln (u), K 2n (u )), ai < u < a 2 } , n = 1,2,... 

be a sequence of Markov processes with values in R 2 whose trajectories belong with probability 
1 to the space D[ ai a 2 ](R 2 ) of cadlag functions on [ 01 , 02 ]. 

Theorem 9. If the finite-dimensional distributions of {K n (u),ai < u < a 2 } converge, as n —>• 00 
to the respective finite-dimensional distributions of a process {K(tt),ai <u< 02 } and for any 
£ > 0 

lim lim sup sup P (||K„(«i) - y|| >e| K„(u 0 ) = y) = 0, 

<d-0 n ~^°° 0<Ui— UQ<g ygR 2 

then, as n —>• 00 

£ {K n (u), ai < u < a 2 } ==>- L {K('o), ai < u < 02 } . 

With this theorem in hands, the method of proving the desired tightness is, essentially, the 
same for all three theorems and has been used in various situations in [ 8 ], [T9] and [ 21 ]. For 
instance, to prove the tightness in Theorems [4] and [5] one can apply a combination of argu¬ 
ments used in j 8 ] and [_19|. The main difference with these previous proofs is related with the 
non finiteness of offspring variances under Conditions ©-Q. Since the remaining parts of the 
proofs follow the reasoning used in [ 8 ] and m we check in this section the tightness for The¬ 
orem [ 6 ] only. In what follows it will be convenient to write P' n ) (£>) for P(£>|Z(n) 7 ^ 0) for any 
admissible event B and to use the notation ||x|| = |xi| + |a? 21 for x = (xi,x 2 ). 


Proof of tightness for Theorem [6] (1) . We will use in this proof ideas from [ 8 ]. As we have 
mentioned, according to Lemma[7]the law P({Z(m, n), 0 < m. < n} E (-)|Z(n) 7 ^ 0 ) specifies, for 
each fixed n an inhomogeneous Markov branching process. We denote its transition probabilities 
by Pl n l(mi, z; m 2 , (■)). In the case under consideration, K„. writes: 


K n (u) = {Z (un, n)|Z(n) / 0} , 0 < u < 1. 
Hence, we will use Theorem [9] for each [0, U] C [0,1). Denote 

C{k) = [z£ : ||z|| < fc} , 

take 0 < uq < u\ < 1 , set 

rrij = Ujn = (1 — e ~ tj ) n, j = 0,1 

and, finally, let 

U(g) = {(ui,u 0 ) € [0,U] x [0, U] : 0 < m - u 0 < g} . 


Lemma 8. Under the conditions of Theorem, \ff{l) for any fixed k and 0 < U < 1 

lim lim sup sup P^ n ^(Z(mi;n) / z|Z(mo; n) = z) = 0 . 

9±° n ~ > ’ 00 U{g)z£C(k) 
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Proof. Let be fixed in N. Notice that when a jump occurs in the process {Z(m, n),0 < 
m < n}, the process |Z(.,n)| + Z±(.,n) is incremented by one. Hence if mo < m\ < m 2 , then 
Z(mo,n) 7 ^ Z(mi,n) implies Z(mo,n) 7 ^ Z(m 2 ,n). Adding Markov property we get 

sup sup P^(Z(mi;n) / z|Z(mo;n) = z) 

U(g) zGC(fc) 

< sup sup P^ n - l (Z(mo + gw, n) / z|Z(mo; n) = z) 

w(s) zeC(fc) 

= sup sup p( ra - m °)(z(gin; n — mo) 7 ^ z|Z(0; n — mo) = z). 
uog[o,i/] z ec(fe) 

According to Lemma [7] and (O 


1 — lim E 


Zi(gn,n—mo) i 
3 i 


Z(n — mo) / 0 


= lim 

n—>• oo 


Qi(gfn; 1 - (1 - si)Qi(n - m 0 - gn)) 


Qi(n - m 0 ) 


= (1 -uo) 1,ai (g+ (1 - si) ai (l-u 0 -g)^ 


— l/o-l 


In particular, 


lim P (n mo) ( 0 ,ei; 5 n,{ei}) = 1 -g/(l -u 0 ). 

n—>-oo 

Therefore, for any e > 0 and (uo,ui) G U(g) there exists no such that for all n > no 

lim P (n) (m 0 ,ei;mi, {ei}) > 1 - g/( 1 - u 0 ) - £ > 1 - g/( 1 - U) - s = 1 - C 15 - e, (58) 


and 


lim lim sup p(”'(mo, ei; mi, {ei}) = 1 . 

gl On^°c u(g) 

Similarly, by Lemma [7] and Theorem 02(1) with a = g/(l — uq) we get 


(59) 


lim E 

n—>-oo 


,Z( mo +9n,n)| Z ^ m0)n ) = 


e 2 


= lim Eo 


3 Z (^- m o)| Z ( n _m 0 ) 7^ 0 


= 1 — a 1 / a2 ^> I (1 — s 1 ) y 


a a 


b V 1 — a 


1/“1 


>( 1 -« 2 ) 


1 — a 


l/a 2 ' 


Further, taking into account the boundary conditions (12TT) and that a 1 = 1 + a 2 yields 


lima 1 /“ 2 '0 ( (1 — s\)— ( 
<40 \ b \ 


a a 


i/«i 


Besides, for 


we have 


At = 


a a 


1 — a 

l/ai 


b \1 — a 


, (1 — s 2 ) 


and At = 


1 — a 


1 — a 


1 / 012 s 


l/a 2 


= 1 - S 2 - 


n —^00 

Note that in view of (12111 


lim p(n-mo) (0)e2 . 5n){e2}) = (i_ a )-l/«2 ^’(^ 2 ) 

oX 2 


(Ai,A 2 )=(ApA^) ' 


dip (Ai, A 2 ) 
d\ 2 


(Ai,A 2 )=(ApA*) ^ 1 C 2(g), 


where C 2 (g) —>• 0 as g | 0. This implies, in particular, that for any e > 0 there exists no such 
that, for all n > no 

P (n) (m 0 ,e 2 ;mi,{e 2 }) > (1 - a ) _1/ “ 2 (1 - C 2 (g)) - e = 1 - C 3 (g) - e 


(60) 
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where C 3 (g) —> 0 as g | 0. In particular, 

lim lim sup Pl”'(mo, e 2 ; mi, {e^}) = 1 . (61) 

gl 0n->oo y(j) 

Using (T59li - (I 6 TT) . the branching property, the decomposability of the process and the positivity 
of the offspring number of each particle of the reduced process, we have for all mi = u\n > 
m o = uqu and z E C(k), 


lim lim inf inf P' n Tmo, z; mi, izj) 
gio u(g) 


2 

lim liminf inf TT(P^(mo, e,-; mi, je,-}))^' 
nitno n—too U(g) 1 \ J J 

J=1 


2 

> limliminf inf TT(P^(mo, e,-; mi, {e,}))^ = 1 . 
gio n-KJO U(g) L1 , J J 

3 =1 


This implies the claim of the lemma. ■ 


Lemma 9. Let the conditions of Theorem^ 1) he valid. If mj = Ujn, j = 0,1,2 and 0 < uq < 
u\ < U 2 < 1 with ui—uo< g, then there exists C^{g) —>• 0 as g f 0 such that for any s > 0 there 
exists no such that for all n > no and z GC(k), 


P^(Z(nii, n) = z 


Z(mo,n) = z; ||Z(m 2 ,n)|| < k) > 


(mo, z; mi,{z})pW (mi, z; m 2 ,C(k)) 
P( n )(m 1 ,z;m 2 ,C(fc))+ (C^g) + 2s) k 


Proof. We have 

P^(Z(mi,n) = z |Z(mo,n) = z, ||Z(m 2 ,n)|| < k) 


P (n) (m 0 ,z;mi,{z}) 


P^(jni,z; m 2 ,C(k)) 
P( Tl )(m 0 , z; m 2 ,C(k))' 


In view of (15Kl) - (l60l) . there exists no such that for n > no, 

P (n) (m 0 ,z;m. 2 ,C(fc)) = ^P (n) (m 0 ,z;m 1; {z'})P (n) (mi, z'; m 2 ,C(A;)) 

z' 


< 1 - P (, ^(m 0 ,z;mi,{z}) + P^(m 1 ,z;m 2 ,C(fe)) 

< 1 - (1 - C 3 (g) - s) k (1 - Cig - s) k + P (n) (m l5 z; m 2 , C(k)) 

< (Cip + C 3 (c?) + 2e) A; + P (n) (m 1 ,z;m 2 ,C(fc)). 


Hence the needed statement follows. ■ 


Let U E (0,1) and 5 E (0, U ) be fixed and, in addition, [7 + 5 < 1. 

Lemma 10. Under the conditions of Theorem. [H( 1) for any fixed k and mo = non, mi = 
uin, m 2 = (U + 5) n we have 

lim lim sup sup P^(Z(mi;n) / z | Z(mo;n) = z, ||Z(m 2 ,n)|| < k) = 0. 

9 l° n ^°°W( 9 ) zeC(fc) 

Proof. Let us first notice that from (1551) we can deduce that if w = (w±,w 2 ) < z = (zi,z 2 ) 
(where the inequality is understood componentwise) then, 


P W (m 0 ,w;mi,{w}) = 




Wl ( j(m 0 ,m i,n) 


3si 


| si=0 


( s ) 




<9s 2 


s =0 


> 


■J 




(»i) 


2:1 / T(mo,r77.i,n) 

/ «^9 


9si 


| 51=0 


( S ) 


dso 


s =0 ] =P W (m 0 ,z;mi,{z}). (62) 
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By Lemma El Equations (l58l) and (l60T) . for m o = non, m\ = u\n and z € C(k) 


p(")(Z(mi, n) = z|Z(mo, n) = z; ||Z(m 2 , ra)|| < k) 


>(1 -C 2 (g) 


e) k (1 - Cig 


P^(mi,z; m 2 ,C(k)) 
P( n )(m 1 , z; m 2 ,C(fc))+ ( 64 ( 5 ) + 2e) k 


Using the decomposability hypothesis, the Markov property of the reduced processes, the in¬ 
equality — mi < (U + 5) n and (1621) we obtain 


2 

P (n) (m 1 ,z;m 2 ,C(fc)) > P ( ' ,) (m 1 ,z;m 2 ,{z}) = JJ(P (n) (m!, e^; m 2 , {ej}))^ 

j=i 

2 2 k 

> U(P^(m 1 ,e J ;m 2 ,{e J })) k > (pW(0,e,; (U + 6) n, {e,})) . 

3 = 1 j =1 

It follows from Subsection 18.31 that 

2 k 2 

limn ( p(rx) (°, ey; (t/ + 5) ri, { e y>)) = J] P fc (W(-log(l-l/-<5)) = ej | W(0) = e^) = H > 0. 

3 = 1 i =1 

Hence we get 


lim inf inf P n (Z(mi; n) = z|Z(mo; n) = z, ||Z(m 2 , n)|| < k) 

n^-ooU(g) zgC(fc) 

> (1 - C 2 (g) - s) k (1 - C ig 


e) k 


B 

B-\- ( 04 (g) + 2c) k 


Letting first g | 0 and then e | 0 completes the proof of the lemma. ■ 

Combining the statement of Lemma [TO] and Theorem El about tightness and the finite dimen¬ 
sional convergence established in Section T8.3I point (1), we obtain the following statement: 


Corollary 1. Under the conditions of Lemma\ 10 \ 


C 


jz(rm, n), 


0 < u < U 


|Z(m 2 ,n)|| < k, Z (n) / 0 j 

=k £ (0jl ) {W (u),0 < u < U\ ||W(C/ + (5)|| < k} . 


Final steps in proving Theorem [6](1): Let for 0 < U < U\ = U + 6 < 1 

P(")([/;(•)) = pW({Z(im,n),fl < u< U} € (•)), 

p( n > k )(U, (■)) = P( n )({Z(un,n),0 < u< U} G (-)| \\Z(Um,n)\\ < k), 

p( n > k )(U, Ui] (■)) = pW({Z(un,n),0 < u< U} € (-)j \\Z(Um,n)\\ > k) 

V(U-(-)) = P({W(u),0<u<U}€ (•)), 

V^(U,U i; (-)) = P({W(u),0 <u<U} e (-)| ||W(C/ 1 )|| < A;). 

Then, for a continuous real function -0 on D[ 0 ) p](Z+) such that 1^1 < q for a positive q we have 
J i/;(x)P {n) (U-dx) = p( n) (||Z([/m,n)|| > k) J ^(x)P^ k \U,U r ,dx) 

+P^(\\Z(U l n,n)\\<k) J 'ip(x)P^’ k \U,U 1 -dx). 

For the first summand we get 


lim sup P <n )(||Z(L r in,n)|| 

n—>00 


>k) j ^(x)P^’ k \u,Ur,dx) 

< glim sup P^(||Z([/in,n)|| > k) 

n—>00 


9P(I|W(17 1 )|| > k) 


o(l) 
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as k —>• oo by the properties of W(-). On the other hand, letting first n —> oo and then k —>• oo 
we obtain from Lemma [TO] and Theorem [9] 


lim lim P n (||Z(C/in,n)|| < k) / ^(x)P (fc) (U,Ui; dx) 

fe—K>o n—>-oo / 


= lim P (0 < ||W(C/i)|| < fc) [ i/>(x)V (k) (U,Ui;dx) 

k—too J 

ip(x)V(U,U 1 -,dx)= [ ip(x)V(U-,dx). 


= lim / 

k ~>°° -v{0<||W(C/i)||<fc} 


Thus, 


lim / ^(x)P^(Z((-)n, n) € dx) = / i/j(x)V(U;dx) 

n^ooj J 

for any bounded continuous function on D[ 0 t[ /](Z^_) proving Theorem [U[l). 


Proof of tightness for Theorem [6] (2). For 0 < t < oo and 1 <C h{ri) <C n let 

Z(„;f) = (7M),+«;*)) = ■ Z i(" ~ Z2( " ~ tt( " > ’" ) ) ' 

For 0 < ti < to w e write 

Zi(n — t\h(n),n) = Z^(n — t\h{n),n-,to ) + zj 2 ^(n — tih(n),n;to) 

where Z^(n — tih(n),n; to), * = 1, 2 is the number of type 1 particles in the reduced process at 
moment n — t\h(n) with type i ancestors at time n — toh(n). Set 


— (+il,+i 2 ) 

and, for e > 0 


( Q2i(n) Q2(n) \ 

\nQi(toh(n)) Zl ’ Q 2 (t 0 h(n)) / ’ 


T ni 


Qi(tih{n)) 
Qi(t 0 h(n)) ’ 


i G {1,2} 


Ani(t) = £nQi(th(n))/Q 2 i(n), A n2 (t) = eQ 2 (th(n))/Q 2 (n). 


Note that under our conditions A n j(t) A oo as n -> oo. With these notations we have 


| Z(n; t\) — z n > 3e| C | \Z 2 (n - tih(n),n) - r n2 z 2 \ > A n2 ( 

U | Z^\n - t\h(n),n\ t 0 ) - r„izi > A„i(ti)| U |zj 2) (n - t\h{n), n; t 0 ) > A n i(ti)| 
To go further we need the representation 

Q 2 (n) 


Z 2 (n;ti) = 


Q 2 {t 1 h(n)) 


Z2(n—t 0 h(n),n) 

^ ^ Cfc(+*o,ii), 

fc=i 


where Cfc( n !^o,^i) = C( n Uo,U) and C( n > to Ti) is the number of type 2 particles in the reduced 
process at time n — t\h(n) being descendants of a particle of the reduced process existing at 
time n — toh(n). Observe that £(n;to,ti) A Z 2 (n — t\h(n)) and EZJ(n — t\h{n)) < oo for any 
7€ [1,1 + a 2 ) in view of Condition ([2]) and Theorem 2.3.7 in [ 12 ] . Notice that 

E [C(n;to, ti)\Z 2 (n - t 0 h(n),n ) = 1] = E 2 [Z 2 ((t 0 - h) h(n),t 0 h(n))\Z 2 (t 0 h(n)) > 0] 

= Q2(hh(n))/Q 2 (t 0 h(n)) 

and that {Cfc(n; to, t\), k = 1, 2,..., Z 2 (n — toh(n),n)} are iid given Z 2 (n — toh(n),n). Using von 
Bahr-Esseen inequality m Theorem 2) for 7 € (1,1 + a 2 ), Markov branching property and 
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Lemma O we obtain that 

(A n 2 (4 )) 7 P (n) ( | Z 2 (n - tih(n),n) - r n2 z 2 \ > A n 2 (4)|Z(n; t 0 ) = s 

| Z 2 (n - tih(n),n ) - r n 2 z 2 | 7 |Z(n;t 0 ) = z n ,Z(rc) / 0 


< E 


< 2E 


22 


y: |Cfc(ra;t 0 , 4 ) - ^n2| 7 |Z(n;t 0 ) = zn, Z(n) 7^ 0 


Lfc=l 


= 2 ^ 2 E 2 [|Z 2 ((t 0 - 4 ) h(n),t 0 h(n )) - r n2 | 7 |Z(t 0 /i(ra)) / 0 ] 
/ 2 7 


< 2^ 2 


E 2 \Z 2 ((to - 4) /z(n), t 0 h(n))} + 2 7 r 7 2 j , 

(63) 


VQ2l(to^(™))' 

where in the last inequality we have applied, for a, b E R 2 

|a — 6| 7 < max(| 2 a| 7 , |26| 7 ) < 2 7 (|a| + | 6 |). 

Further, let 

4 (t 0 , 4 ;n), k = 1 , 2 , ...,Z 2 ((t 0 - 4) /i(n)) 

be the indicator of the event that the fcth particle belonging to generation (to — ti) h(n) has 

descendants in generation toh(n ) and Ik (4,4; n) = / (to, 4; n). Now, again by (l63li and the 
Bahr-Esseen inequality 


2 7 E 2 [|Z 2 ((t 0 - 4 ) /i(n),t 0 ti(n))| 7 ] 

Z2((to—ti)h(n)) 

= 2 _ 7 E 2 ^ (4 (t 0 , 4 ;n) - Q 2 (4/i(n))) + Z 2 ((t 0 - 4 ) h(n))Q 2 (tih(n)) 

k =1 

Z 2 ((t 0 -t 1 )h(n)) 

< E 2 (4 (to,4;ra) - Q 2 (tih(n))) 

k =1 


71 


71 


+ E 2 [Z 7 ((t 0 - 4) h(n))Ql(tih(n))] 


< E 2 [Z 2 ((t 0 - 4) h(n))} E 2 [|l{z 2 (tih(n))>o} “ <32(4ti(n))| 7 ] 

+ Q](4/t(n))E 2 [ZJ((t 0 — 4) h(n))] . 

Observe now that 

E [| 1 {^ 2 (*ift.(ri))>0} - Q 2 {tih(n ))\ 7 ] 

= Q 2 (ti/i(n)) |1 — Q 2 (ti^(n))| 7 + (1 - Q 2 (4/i(n))) Ql(tih(n)) < CQ 2 (tih(n )) 
and by Lemma 11 in [18] 

E 2 [Z 7 ((t 0 - 4) M n ))l < CQ^ 7 ((to - 4) M n )) < CQ 2 -7 (t 0 h(n)). 

As a result we get 

E 2 [ZJ((t 0 - 4) h(ri),t 0 h(n))] < CQ 2 (tih(n)). 

This shows that, for each fixed pair 0 < p\ < p 2 < 00 


P (n) ( |Z 2 (n - tih(n),n) - r n2 z 2 \ > A n 2 ( 4 )|Z(n;t 0 ) = z n ) < C4 2 ( — ) < C ( 

which goes to 0 as n 00 uniformly in (to, 4 ) € \pi,p 2 ] and 

Z 2 <CK n2 (p 1 ) or z 2 = z 2 /A n2 (t 1 ) < Ci. 


1 \7-l 


A n (pi) 
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Further, for 7 E (1,1 + 07 ) we have 


p(n) 



t\h{n),n\ to) 


r n izi 


> Ani(ti)| Z(n; t 0 ) 



< 


1 


We now use the representation 

hi) 


E 2 


(1) 7 - 

Z{ (n — fih(n),n;f 0 ) - r n iz± | Z(n;t 0 ) = z„,Z(n) / 0 


Z\ (n—toh(n),n) 


Z^\n - tih(n),n-,t 0 ) = y r]k(n]t 0 ,t 1 ) 


k =1 


where 7 ifc(n; to, H) = 7 (n;to,ti) and 7 (n;to,ti) is the number of type 1 particles in the reduced 
process at time n — t\h{n ) being descendants of a type 1 particle existing in the reduced process 
at time n — toh(n). Observe that 

E [ 7 7 (n; to, ti)] < Ei [ Z'l (n — t\h{n))\ < 00 

for any 7 € [1,1 + cei) in view of Condition ([2]) and Theorem 2.3.7 in [12]. Similarly to the 
previous estimates 


E 


Z\ (n—toh(n) ,n) 

y \Vk(n; to, t\) - r„i | 7 |Z(n;t 0 ) = z n , Z(n) / 0 

k=1 


y 2z\ Ei 

— 2zi~Ei 


V(n-,t 0 ,ti) - 


Qi(tih(n)) 


Qi{toh{n)) 

Zi{{to - ti) h(n),t 0 h(n)) - 
Hence, for each fixed pair 0 < p\ < p 2 < 00 we get 


\Zi(n — toh{n),n) = 1 

\Zi(t 0 h(n)) > 0 


Qi{t\h{n)) 7 


Qi(t 0 h(n)) 


< C Zl . 


p( n ) ( Z^\n - tih(n),n;to) - r n izi > A„i(ii)|Z(n; t 0 ) = z) < 


1 


A n l(tl) 


Z\ —>• 0 


as n —>• 00 uniformly in (to,ti) € [pi,P 2 ] and z \ + Ch? nl (jpi) or £1 = zi/A^^i) < Cj. 
Finally, using the branching property, we get 


A n i(fi)P (n) (zj 2) (ra - tih(n),n;t 0 ) > A n i(ti)| Z(n;t 0 ) = 


< E 


Z{ 2) (n - tih(n),rr, £ 0 )| Z(n;t 0 ) = z n , Z(n) / 0 
= z 2 E 2 [Zi((t 0 - h) h(n);t 0 h(n))\ Z (t 0 h(n)) / 0 ] 

A 2 i (t 0 - h) h{n)Q\{t\h{n)) 


= Z 2 - 


Q 2 i(t 0 h(n)) 


< C 


Z2 


Anl(tl) 


Thus, for each fixed pair 0 < p\ < p 2 < 00 

p(n) 7 Z (2) 


(n - t 1 h(n),n;t 0 ) > A n i(ti)|Z(ra; t 0 ) = z n ) < Cvr 

7 A nl 


- 2-2 


(tl) 


as n -A 00 uniformly in (to, ti) € [pi,p 2 ] and z 2 = 2 2 /A re i(ti) < (7]. 

The remaining part of the proof follows the line of proving Theorem ( 6 j[ 1) and is based on the 
one-dimensional convergence established in Theorem 0(2). We omit the details. 
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10. Most recent common ancestor 


We present in this section the derivation of the laws of the MRCA. They are direct conse¬ 
quences of Theorems |4]|6j 

Point (1): It follows from Theorem [4] that, given Condition <[ 8 |) particles of type 2 are always 
present in the limiting process. Thus, the MRCA should be of type 2. Since, given Z (n) 7 I 0 
the event = 1 } has a nonnegligible probability only if lim^^oo ran -1 < 1 , to find the 

distribution of the MRCA in this situation it is sufficient to evaluate the derivative with respect 
to s 2 of the limiting probability generating function (fT4D at point S 2 = 0. This gives, as desired 


lim P2 ( Z 2 ( an,n) = l|Z(?r) 7 ^ 0) = lim P2 (/ 3 n > em|Z(n) 7 ^ 0) = 1 — a. 

n—>• 00 n—> 00 


Point (2): According to the offspring generating function of the limit process, given in (|24j) . 
the law of the offspring vector (£ 21 ; £ 22 ) produced by a type 2 individual satisfies: 

Qo 1 

p( 6 i = 0,62 > 2) = —- and P(&i = 1, £22 = 0) = —-. (64) 

1 + «2 1 + «2 

Hence there are two possibilities: 

• Either the initial type 2 individual generates several type 2 individuals at the firth birth 
event. Then the time of the branching corresponds to the death of the MRCA. 

• Or the initial type 2 individual gives birth to a type 1 individual at the first birth event. 
Then the process does not evolve anymore as the type 1 individuals are immortal and 
sterile. We will see that in this case the MRCA birth time is of order larger than g*{n). 

From ([MD we get that 

P(Yi(i) + Y'lit) = 1 for any t > 0) = --. 

1 + «2 

Moreover, from points (2)-(3) of Theorem [2J we see that for g*{n ) <C m <C n 

lim P 2 (Zi (m,n) = l|Z(n) / 0) = —-—, 

n—>oo 1 + CI 2 


and for a E (0,1) 

lim T* 2 (Zi( an, n) = l|Z(n) / 0) = --. 

n-t 00 1 + Q 2 

As the process Zi(m,n) is non-decreasing with m, we deduce that for g*(n) m ^ n and 
a € ( 0 , 1 ), 

lim P 2 (m <C /3 n < an|Z(n) / 0) = -—-—. 

n-H3 o 1 + «2 

If, on the contrary the initial type 2 individual gives birth to several type 2 individuals, which 
happens with a probability 02/(1 + 012 ), the time of this first branching will be an exponential 
variable with parameter (1 + 0 . 2 )! 0 . 2 - Converting this distribution in terms of the limit law for 
the MRCA yields 


lim P 2 (/?n < tg*(n)\Z(n) / 0 ) = 2 (1 - e t ( 1 + a ^)/ a ^). 

n-too 1 + 02 

This ends the proof of Theorem 0(2) . 

Point (3): To find the distribution of the MRCA birth time under Condition (fTOl) we use 
the offspring generating functions of the process |W(t) ,t> 0}, which are given in (12711 and 
(1281) . Equation (l28j) implies the following properties for the law of the offspring vector (£ 21 ,^ 22 ) 
produced by a type 2 individual of W(-): 


IM&I = 1,^22 = 0) 


(702^21 

6((1 + a 2 )^“ 2 


■jy) P 2 ( 6 i = 0,^22 > 2 ) 


«2 b a2 

(1 + 02 ) 6“ 2 — 1 


(65) 
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Moreover, we see from the definition of W±(-) and (1271) that a type 1 individual dies after an 
exponentially distributed time with parameter 1 and gives birth to two type 1 individuals at 
least. Hence if we denote by T the time of the first branching there are two possibilities: 

• Either the initial type 2 individual gives birth to several type 2 individuals. In this case 
T will be an exponential variable with parameter k, where n has been defined in (1261) . 
According to (1551) this happens with probability 

0:2 6" 2 b a2 

(1 + 0:2)6" 2 — 1 K 

• Or the initial type 2 individual gives birth to one type 1 individual. Then this type 
1 individual has an exponential distributed life-time with parameter 1. In this case, 
T will be the sum of an exponential variable with parameter n and of an independent 
exponential variable with parameter 1. According to (1651) this happens with probability 

(JO2A12 _ 0A21 

6((1 + 02 ) 6" 2 — 1 ) 6 k 


Hence we get 


P(T < t) 



K 



Q~Ai 2 

bn 



1 


K — 1 


e 


— tit 



By using several times the equality 6 1+ " 2 — 6 = craoA 12 given in (1191) . we finally obtain 

1 

—p. L — 

l + 0 2 


P (T <t) = 1 - 


_g —t _ t((l+a2)f>“ 2 -l)/Q2 


1 + o 2 

Then by applying Theorem [6]) 1), we get for every a G (0,1): 


lim P {p n < an) = P(T < — ln(l — a)). 

n—>00 

This ends the proof of Theorem 0)3). 
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